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Abstract 

. In 1976, Shlyakhter showed that the Sm data from Oklo results in 

the upper bound on the time- variability of the fine-structure constant: 
\a/a\ <,10 _17 y _1 , which has ever been the most stringent bound. Since 
the details have never been published, however, we recently re-analyzed 
^ \ the latest data according to Shlyakhter's recipe. We nearly re-confirmed 

his result. To be more precise, however, the Sm data gives either an upper- 
Q\ ' bound or an "evidence" for a changing a: a/a = — (0.44±0.04) x 10 _16 y _1 . 

. A remark is made to a similar re-analysis due to Damour and Dyson. We 

also compare our result with a recent "evidence" due to Webb et al, ob- 
tained from distant QSO's. We point out a possible connection between 
this time-dependence and the behavior of a scalar field supposed to be 
responsible for the acceleration of the universe, also revealed recently. 



o 
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1 How did constants cease to be constant? 

o 

More than 60 years ago, Dirac startled the science community by saying that the 
gravitational constant G is not a true constant, but decays like the inverse of the 
cosmic time t: 1 

G(t)~r\ (1) 



from which follows that the rate of the fractional change G/G should be like 

§ ~ -r 1 ~ -io-V\ (2) 

because the present age of the universe is believed to be somewhere around 10 10 y. 

I do not know how many physicists were really convinced with Dirac's argu- 
ment, based on the "Large- numbers Hypothesis." No matter what philosophy he 
had started with, however, it is clear that he opened up a new view of looking at 
the very fundamental aspect of laws of Nature. Some people tried hard to prove 
or disprove his result. But the number is so small that almost any noise would 
easily wipe out the real effect if there is any. 
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So far no solid evidence for the changing G has been reported. We have heard 
only about upper bounds, like 



These are already below Dirac's prediction. The search for the changing G may 
appear to have lost the motivation. We have now, however, "unified theories," 
attempts of unifying physics of gravitation and of microscopic world. According to 
the theoretical models, including string theory as one of the best known examples, 
time-dependence of various coupling constants should be expected generically for 
the following reasons: 

• Truly constant coupling constants in the theories at the deeper level in 
higher-dimensional spacetime. 

• At another level we have a realistic theory in 4 dimensions, where the ef- 
fective coupling constants to be measured are derived from those in the 
fundamental theories through the process of compactification. 

• In this process, some of the scalar fields, notably a dilaton, comes into play. 
The scalar fields are expected to change slowly as the universe evolves, and 
so do the effective coupling constants. 

In this context, the way the coupling constants change may be different from 
the way Dirac envisioned. Obviously coupling constants are not only confined 
to the gravitational coupling constant, but also include other coupling constants, 
notably the fine-structure constant a = e 2 /(47ihc), the square of the electromag- 
netic coupling constant and its strong interaction analogue a s . The result of the 
past investigations on their time variability is summarized: 

Most of them are upper bounds, except for the last entry on the recent result 
due to Webb et al, who claim the first "evidence" of the changing a. 2 As another 
point to be noticed in this Table, we find that the "Oklo phenomenon" gives an 
upper bound several orders of magnitude better than any other results. 

2 Oklo phenomenon 

Oklo is the name of a uranium mine in Gabon, West Africa, but should be 
remembered as a place where the evidence of "natural reactors" was discovered 
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Table 1: Results for time- variability of a and a s . 



and verified for the first time in 1974. 3 By natural reactors, we mean that self- 
sustained nuclear fission reactions occurred naturally about 2 billion years ago, 
lasting some million years. This was a rediscovery of the theoretical prediction 
due to Kuroda in 1955. 4 

Though it might sound rather strange, this phenomenon took place for several 
reasons. The most crucial among them is the fact that the abundance of 235 U 
must have been ~ 3%, much higher than the present value 0.7207% due to the 
difference between lifetimes of 235 U and 238 U. Since 1974, extensive studies have 
been made on the remnants of fission products in the so-called reactor zones in 
Oklo. 

In 1976, Shlyakhter pointed out that observation of the isotope 149 Sm could 
be useful in determining how much different the nuclear reactions 2 billion years 
ago were from what they are today. 5 The natural abundance of this isotope is 
13.8%, while the observation in Oklo showed much smaller values. This deficit 
can be understood in terms of the consumption due to the neutron absorption 
reaction: 

n+ 149 Sm-^ 150 Sm + 7 , (3) 

that had occurred inside natural reactors. 

This scenario was verified by using the data obtained from today's laborato- 
ries. This implied that the parameters of the reaction (|3p should not have been 
very much different from what they are today. This conclusion was made even 
more precise by noting that this reaction is dominated by a resonance lying as 



3 



low as E r = 97.3meV, which is smaller by more than 7 orders of magnitude than 
a typical energy scale ~ lMeV in most of the nuclear reactions. The half-width 
is also very narrow; 61.0meV. 

In view of the smallness of E r , we should consider that the presence of this 
resonance is a consequence of a nearly complete cancellation between two in- 
teractions; repulsive Coulomb energy and the nuclear force due to the strong 
interaction, leaving a tiny residual effect. Suppose we fictitiously change the 
value of a. No matter how slight the change Aa might be, the residual change 
AE r of the resonance energy could be substantial, particularly in terms of the 
relative change AE r /E r . The same must be true then for the cross section <t 149 
itself. 
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Figure 1: Effective cross section for n + 149 Sm -^ 15 ° Sm + 7, for thermally averaged neutrons, 
plotted against AE T . 

This is an excellent example of an enhancement mechanism, only by which 
one can measure something very small. In Fig. 1, we plotted the cross section 
<3"i49 as a function of AE r , also assuming that the neutron flux is in the thermal 
equilibrium of the temperature T. The hat attached above a implies that the 
Boltzmann distribution is normalized to a particular temperature T = 20.4°C, 
as a convention in most of the Oklo analyses, hence the ordinary cross section 
a is given by cr(T) = ^ (n / 4)(Tq/T)ct(T) . We see how sensitively a depends on 
AE r . If E r 2 billion years ago were to have been smaller than the present value 
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by as little as lOmeV (AE r = — lOmeV) and T = 20°C, for example, the cross 
section would have been larger by about 10%. 

The trouble is, however, his paper in 1976 was too short to give any details 
on what the numerical result was like as well as of what quality his data was. It 
is rather difficult to know how reliable his final result was. We finally decided 
to re-do what he did previously. We reconstructed his formulation based on his 
unpublished paper. 6 

We soon realized that a major source of error comes from the following 
situation. 7 Each of the natural reactors must have come to an end in millions 
of years. Even after this termination, some amount of the isotope might have 
migrated from outside into the core part. This amount has nothing to do with 
the nuclear interaction we are interested in. We tried to minimize this contami- 
nation due to weathering and related phenomena as much as we could. For this 
purpose we chose five samples taken from deep underground with a great care of 
geologist's expertise. In this way we obtained the data of quite good quality. 

We analyzed the evolution equations for the system of 235 U, 147 Sm, 148 Sm, 149 Sm, 
obtaining 

5"i49 = 91 ± 6kb. 
We used several ways to estimate the temperature, finding 

T = 200 - 400°C. 

Applying these ranges to Fig. 1, we found two separate ranges:[|] 

!9 ± 11 meV, for the right-branch range, 

—97 ± 8 meV, for the left-branch range. 

Two ranges were obtained because a horizontal line corresponding to the ob- 
served cross section intersects the theoretical curve of the cross section for each 
temperature at two points, at both sides of the peak. 

The right-branch range covers AE r = 0, giving a null result with an upper 
bound as usual, whereas the left-branch range indicates, as it is, that the value 
of AE r was different from today's value by 12 standard deviations. 

We use a theoretical analysis: 8 

Aa _ AE r 

^The result for the right-branch range replaces Eq. (27) in reference 7, where the erroneous 
lower end was shown inadvertently. 
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where the mass scale M c is estimated to be 

M c « -l.lMeV. 
This allows us to translate AE r to the change of a: 

Aa _{ -(0.8±1.0)xl0- 8 , 
a ' \ (0.88 ±0.07) x lO" 7 . 

By dividing by —2 x 10 9 y we finally arrive at 

a_ | (0.4±0.5)xlO- 1 V 1 , 
ol [ -(0.44 ±0.04) x KT^y -1 - 

As it turns out, the result from the right-branch agrees very well with what 
Shlyakhter obtained in 1976. In this sense we have confirmed his result. The 
agreement to this extent, however, seems rather accidental, because it is unlikely 
that the data as good as ours was available at early years of Oklo investigations. 

We tried hard if we could eliminate the non-null result from the left-branch 
range by using other isotopes. We had results from 155 Gd and 157 Gd for which ex- 
ceptionally low-lying resonances play important roles. Unfortunately, we reached 
the conclusion short of complete, because better "efficiency" of the resonances 
made the residual abundance even scarcer, and left results more "contaminated." 
The null result is certainly favored but not beyond that. 

We also mention, on the other hand, that the same kind of re-analysis had been 
made by Damour and Dyson, who used, however, older and more contaminated 
data. 8 As a result, they reached a much wider range of AE r , essentially covering 
our two ranges together without separation between them. 

3 Result from distant QSO's 

By a recent analysis of the spectral lines from distant QSO's, Webb et at. obtained 
the first evidence for the changing a? They exploited relativistic corrections to 
the spectrum having another a 2 -dependence compared with the Balmer series 
spectrum. In this way they enhanced the number of lines to be analyzed by a 
factor 10, hence the statistics by 100, compared with the measurements of alkali 
doublet splitting made so far. By using the data from 49 QSO's with the look- 
back time ranging from 40% to 90% of the age of the universe ~ 1.4 x 10 10 y, they 
concluded 

Arc 

— = (-7.5 ±1.8) x 10- 6 , 
a 
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indicating a non-zero result beyond 4.3a. 

Our Oklo result is at least 2 orders smaller, but the look-back time is also 
smaller; xs 15% in the above sense. Combining them all it looks as if the time- 
dependence as a function of the cosmic time was nonmonotonic, or even oscillat- 
ing. 

4 An accelerating universe 

Recent analyses of Type la supernovae finally concluded that the universe is 
accelerating, 9 described best in terms of a small but nonzero cosmological constant 
with the size as given by 

tt A = 0.6-0.7. 

This turns out be in agreement with other kinds of analyses, including the age of 
the universe, large-scale structure, number-counting of distant galaxies and lens 
statistics. We used to have the "first face" of the cosmological constant problem; 
why it is smaller than the theoretically natural value ~ Mp by as much as 120 
orders of magnitude. Now we confront the "second face;" how we can understand 
a nonzero fl\ rather close to unity. This is closely connected with the "coincidence 
problem," why are we so lucky to be in this very special epoch in which Q\ ~ 1? 

A possible reply to the first question is to call for the scenario of a decaying cos- 
mological constant; A(t) ~ t~ 2 ~ 10~ 120 . This has been implemented successfully 
in terms of a version of the scalar-tensor theory. 10 According to this scenario, to- 
day's A is small only because we are old. More specifically, we think it convenient 
to use the reduced Planckian unit system in which c = h = M P (= (8-kG)^ 1 ^ 2 ) = 1. 
The age of the universe t « 1.3 x 10 10 y is given in units of the Planck time as 
logt ~ 60, from which follows p cr ~ t~ 2 ~ 10~ 120 , compared with a theoretical 
estimate A th ~ 1, the easiest way to understand a number ~ 10 120 . 

We studied a prototype Brans-Dicke model with A added to the Lagrangian: 

C = V=9 - ^iTW - A + Matter) > ( 4 ) 

with the scalar field and the BD parameter u, where we use the notations, 
e = Sign(cj) and = (l/4)^~ 1 . The first term on the right-hand side is called a 
nonminimal coupling term which implies that we have an effective gravitational 
constant G e s = (87r£0 2 ) -1 , which is spacetime dependent, in general. 

We emphasize that a conformal transformation ds 2 — > ds 2 = £l(x)~ 2 ds 2 , with 
a local function Q(x), is inherently important in any theory with a nonminimal 
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coupling term. By a conformal transformation we move from one conformal frame 
to another. A special conformal frame in which the Lagrangian is given by (|J) is 
called the J(ordan) frame, while the conformal transformation chosen by 

5V = ^ 2 9*^u, with Q 2 = £(fi 2 , (5) 

takes us to the E(instein) frame, in which the same Lagrangian (Q) is re-expressed 
as 

C = ^Tgl (~R* - ^d^advo - Ae' 4 ^ + L* matte r) , (6) 
where the new scalar field a is defined by 

= r 1/2 e C<7 , C 2 = 6 + eC\ (7) 
with the coupling strength (. Substituting this into the second of (|) yields 

Q = e c<T . (8) 

We notice that the nonminimal coupling term has now be removed and the A 
term acts as a potential, he~^ a . The energy density of a allows an interpretation 
as an effective cosmological constant A e g in the E frame. It has been argued 
frequently that the J frame is a physical conformal frame where we live. Our 
study of the cosmology with A added reveals, however, that the E frame is to be 
preferred with a remedy assumed in the matter Lagrangian. 11 It also turns out 
that the quantum anomaly effect coming from interactions among matter fields 
is crucially important. 

We also derive the (classical) cosmological equation in spatially flat Friedmann 
universe. The attractor solution can be obtained, which yields 

Pa = Acs = \o 2 + Ae~ 4CCT ~ r 2 ~ p m , (9) 

where t* is the cosmic time in the E frame. This equation shows that the effective 
A does decay like t~ 2 , as expected. This is the behavior sometimes called "scaling" 
in the current discussion of "quintessence, 12 " and is a necessary condition for the 
solution of the coincidence problem. 

We find, however, that this behavior is precisely the same way as the ordi- 
nary matter density p* m falls off, as also shown in (§). This situation does not 
correspond to the observational finding which requires, in the present context, 
that p a stays nearly constant at least temporarily, not necessarily forever, thus 
mimicking a cosmological "constant." In order to meet this second face, we make 
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use of a feature called a "hesitation behavior," which is unique to a fast falling 
potential like an exponential potential. This implies that there is a duration of 
time when the scalar field stays nearly constant, thus is expected to act like a 
cosmological constant. We find, however, this "built-in" mechanism is short of 
resulting in such a "large" ~ 0.7. 




20 40 60 

log t. 



Figure 2: An example of the solution of the classical equation for the potential V(a, x). Upper 
diagram: b = In a* (dotted), a (solid) and 2\ (dashed) are plotted against logi*. The present 
epoch corresponds to logi* = 60.1 — 60.2, while the primordial nucleosynthesis must have taken 
place at logi* ~ 45. The parameters are A = l,( = 1.5823, m = 4.75,7 = 0.8, k = 10. The 
initial values at t\ — 10 10 are a\ = 6.7544, a[ — (a prime implies differentiation with respect 
to r = hit*), xi = 0.21, xi = -0.005, pi rad = 3.7352 x 10~ 23 , pi dus t = 4.0 x 10~ 45 . The dashed- 
dotted straight line represents the asymptote of a given by r/ (2(). Notice long plateaus of a and 
X, and their rapid changes during relatively "short" periods. Middle diagram: a* = b' = t*H* 
for an effective exponent in the local power-law expansion a* ~ t°* of the universe. Notable 
leveling-offs can be seen at 0.333, 0.5 and 0.667 corresponding to the epochs dominated by the 
kinetic terms of the scalar fields, the radiation matter and the dust matter, respectively. Lower 
diagram: logp s (solid), the total energy density of the a — x system, and logp* m (dashed), the 
matter energy density. Notice an interlacing pattern of p s and p* m , still obeying ~ t~ 2 as an 
overall behavior. Nearly flat plateaus of p s precede before p s overtakes p* m , hence with Q\ 
passing through 0.5. 

We searched for a remedy on a try-and-error-basis, coming across finally a 
"two-scalar model." We introduced another scalar field, called x, together with 
an interaction potential in the E frame: 

where 

V = A + -m 2 x 2 [1 + 7 sin(«;cr)] , 

with m, 7, and k constants chosen not to be very much different from unity in the 
reduced Planckian unit system. Note that this potential goes into the exponential 
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Figure 4: Magnified view of a (solid) 
and 0.02<I> + 44.25 (dashed) in the upper 
diagram of Fig. |[ Note that the verti- 
cal scale has been expanded by approx- 
imately 330 times as large compared 
with Fig. ||. With the time variable r = 
In , the potential V grows as the mul- 
tiplying factor t\ = e 2r . The potential 
wall for a becomes increasingly steeper, 
thus confining a further to the bottom 
of the potential, as noticed by an oscil- 
latory behavior of a with ever increas- 
ing frequency measured in r. The grow- 
ing V causes x eventually to fall down- 
ward, resulting in the collapse of the 
confining potential wall. The stored en- 
ergy is then released to unleash a. 

potential as before in the limit \ = 0. This potential might look complicated and 
artificial, but turns out to result in many desired features, as one will find in a 
typical example of the solution illustrated in Fig. ^|. See Ref. 13 for more details. 

We fine-tuned parameters and initial values moderately. Around the present 
epoch, log^o ~ 60, the density p s = A eS for the scalar-field system surpassed 
the ordinary density p vm as shown in Fig. |, giving Qj^ = 0.75 with the Hubble 
parameter Hq = 72.9km/sec/Mpc. The coincidence problem is now less severe 
because a "mini-inflation" is one of the repeated events, instead of a once-and- 
for-all occurrence in the whole history of the universe. 

Figs. 3,4 are magnified view around the present epoch. In particular, Fig. 4 
shows an oscillatory behavior of cr, which might be related to the time variation 
of a, as will be discussed. 



59 59.5 60 60.5 61 

log t. 

Figure 3: Magnified view of \ogp s 
(solid) and logp* m (dashed) in the 
lower diagram of Fig. || around the 
present epoch. Note that p s is very flat 
in this diagram extending back to the 
past of z — 5. 2 — 6.9 for the assumed 
age (1.1 - 1.4) x 10 10 y- 
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5 Time-dependent fine-structure constant? 



We consider QED as matter part of the Lagrangian (Q) in the J frame, by intro- 
ducing a charged spinor field ip and the Maxwell field: 

^matter = V=9 {$ ' + ™4 + ™) 4 ~ ^) , (10) 

where the charge e is assumed truly constant in this J frame. 

We apply the conformal transformation. In the E frame we include the loop 
integrals which are divergent, hence requiring regularization. For this purpose we 
assume .D-dimensional spacetime. We stay in D dimensions until we go to the 
limit D — > 4 at the end of the calculation. 

To be noticed is the fact that the Maxwell theory is conformally invariant in 
4 dimensions. It then follows that the electromagnetic field in D dimensions 
is transformed into A^^ defined by 

A*^ = Vl A^, (11) 

where d = D/2. It also follows 

e, = Q d ~ 2 e = ee {d - 2Ka w e [1 + (d - 2)(a] . (12) 

It appears that e remains unchanged for d — *■ 2. According to the spirit of 
dimensional regularization, however, we leave d — 2 as it is at the moment. (|i~2j) 
implies that a change of the field a induces a classical change of the electric charge 

Ae* = C(d - 2)e*A(7. (13) 

We consider the electromagnetic vertex of ip, inserting a if) loop in the photon 
self-energy part, which develops a pole (2 — This cancels the zero (d — 2) 

in ([13|) , leaving a finite contribution, in the same manner as we obtain quantum 
anomalies. In this way we arrive at 

^1 = ^lAa w 1.2 x 10- 3 Aa. (14) 
a* 2ir 

From the example shown in Fig. ^, we find Acr w 5 x 10 -2 during the time 
span 10 59 — 10 60 corresponding to the look-back time of 90% of the age of the 
universe. Using this on the right-hand side of (|i~4j) , we finally obtain 

^^0.6xl0" 4 , (15) 
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which is found to be an order of magnitude too large compared with the obser- 
vation in Ref. 2. 

We need to select other solutions by changing the initial values slightly, with- 
out changing the result for and H in any appreciable amount. Also we have 
to improve the theoretical calculation with respect to other particles inside the 
loop. 

At this moment, we only say that comparison with the result on Aa/a con- 
strains the cosmological solution much further beyond fitting the accelerating 
universe alone. Details of the analysis will be published elsewhere. 14 | 

6 Concluding remarks 

Based on an assumption in the scalar-tensor theory, we have come to providing 
an explicit implementation of the idea that time-dependence of the fine-structure 
constant ensues from the evolution of the scalar field, which is responsible for the 
accelerating universe. Still even at this preliminary stage, the order-of-magnitude 
agreement between the sizes of the two phenomena seems remarkable. 

As another outcome unique to the scalar-tensor theory, we point out that the 
proposed remedy to the prototype Brans-Dicke model entails the occurrence of 
non-Newtonian gravity featuring a finite force-range and the violation of Weak 
Equivalence Principle probably within the observed constraints obtained so far. 16 
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